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Abstract —In this paper, we propose fusion of dynamic 
TOA (time of arrival) from multiple non-coherent detectors 
like energy detectors operating at sub-Nyquist rate through 
Kalman filtering. We also show that by using multiple of these 
energy detectors, we can achieve the performance of a digital 
matched filter implementation in the AWGN (additive white 
Gaussian noise) setting. We derive analytical expression for 
number of energy detectors needed to achieve the matched 
filter performance. We demonstrate in simulation the validity 
of our analytical approach. Results indicate that number 
of energy detectors needed will be high at low SNRs and 
converge to a constant number as the SNR increases. We also 
study the performance of the strategy proposed using IEEE 
802.15.4a CM1 channel model and show in simulation that two 
sub-Nyquist detectors are sufficient to match the performance 
of digital matched filter. 

Index terms: Time of arrival (TOA), ultra wideband (UWB), 
UWB ranging. 

I. Introduction 

An ultra wideband (UWB) system is based on spreading 
a low power signal in to large bandwidth. There are sev¬ 
eral UWB techniques, impulse radio based UWB (IR-UWB) 
schemes are most popular as they provide better performance 
and complexity trade-offs compared to other UWB schemes 
0-0. Narrow impulse signals that are used in IR-UWB 
schemes yield very fine time resolution and thus, can be used 
for accurate measurement of time of arrival (TOA) m, a. 

Accurate TOA is essential in several applications includ¬ 
ing localization and communication. In localization, when 
the nodes are synchronized, the TOA of the signal can be 
used directly to obtain the range estimate. If nodes are not 
synchronized, TOA estimate is still needed for several ranging 
protocols to estimate the range. Localization information about 
the node can be derived from its range to anchor nodes 0, 
0. In many IR-UWB communication systems, the message 
information is embedded in the location of the IR-UWB pulse. 
For example, pulse position modulation (PPM) and its variant 
modulation schemes. To demodulate the IR-UWB symbols in 
these modulation schemes, TOA estimation techniques can be 
used 0-0. 

IR-UWB system with matched filter can provide optimal 
estimate of TOA for additive white Gaussian noise (AWGN) 
limited channels. The TOA is estimated by finding the peak 
location at which the filter attains the maximum value. 
The impulse response of the matched filter should closely 
approximate the transmit pulse shape. When implemented 
digitally, matched filter requires Nyquist-rate sampling Col, 
no. There are several digital UWB receiver architectures for 


which matched filter is suitable including fl2l - llT4l . Digital 
receiver structure offers several benefits such as flexibility in 
design, reconfigurability and scalability fBl . However, UWB 
signal occupy extremly large bandwidth and thus requires 
high-speed analog-to-digital converters (ADCs). These speeds 
demand the use of interleaved flash ADC ca or a bank of 
polyphase ADCs CD- In addition, the ADC must support a 
large dynamic range to resolve the signal from the strong 
narrowband interferes. These aspects makes digital UWB 
architecture which operates at Nyquist-rate to be costly and 
power hungry |fl4l . IT6l . fT7l . 

Due to the above discussed practical concerns, energy detec¬ 
tors are an interesting alternative as they are less complex and 
can operate at sub-Nyquist rate fT8l - li20il . However, these de¬ 
tectors suffer from the noise due to the square-law device and 
are sub-optimal in AWGN channels. Performance analysis of 
stand alone energy detector for TOA estimation in sub-Nyquist 
rate has been studied in 03, ED, ED- However, how to 
achieve the performance of a digital matched filter using mul¬ 
tiple energy detectors for a dynamic TOA model is to the best 
of the authors knowledge unavailable. In this paper, we first 
show that TOA estimation using digital matched filter achieves 
Cramer-Rao lower bound (CRLB) asymptotically with SNR. 
We show analytically that we can achieve the performance of 
a Nyquist rate matched filter, by using multiple sub-Nyquist 
rate energy detector estimates. We show that number of energy 
detectors needed to meet the matched filters performance is 
high at low SNRs and reduces as SNR increases, and finally 
converges to 4 as SNR increases asymptotically. We propose 
joint fusion and tracking of dynamic TOA using multiple 
energy detector estimates through a suitable Kalman filter 
design. We assess the performance of fused TOA estimate in 
simulation for a dynamic TOA model from multiple energy 
detectors and compare it with the digital matched filter and 
demonstrate the validity of the claims. 

This paper is organized as follows: In Section mi we 
discuss the signal and system model employed. In Section 
Hill we will discuss different types of detectors. Here, we 
will study the energy collection strategy for matched filter 
and energy detectors. In Section IIVI performance analysis 
of these two detectors are studied for a static TOA case. 
Section PvTl discusses the Kalman filter design for joint fusing 
and tracking of the energy detector estimates for a dynamic 
TOA model. Section IVII1 demonstrates the performance in 
simulation. Finally, Section IVIlll details the conclusions from 
the design and results demonstrated. 
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Fig. 1. Joint fusion and tracking structure using N distributed receivers using low cost energy detectors operating at sub-Nyquist rate, tracking a dynamic 
TOA of a moving target. 


II. Model 

A. Signal Model 

The signal model comprises of Nf frames each having a 
unit energy pulse s(i), given by 

N t 

^tr(f) = Y, djs(t - jTf - CjT c ), (1) 

3 = 0 

where each frame is of duration Tf and the frame index is 
represented by j. The chip duration is represented by T c and 
Cj e {0 ... N c } indicates the time-hopping code, dj € {±1} is 
the polarity code, which can be used along with time-hopping 
to smooth the signal spectrum. We consider the frame duration 
longer than the delay spread of the channel. 

A wide range of pulse shapes have been explored for UWB 
applications from rectangular to Gaussian. Gaussian and their 
derivatives, usually called monopulses, are effective due to 
the ease of construction and good resolution in both time 
and frequency. In many cost effective hardware designs, these 
shapes are generated without any dedicated circuits 0, 0, 
1123 1 . (24). For analytical and simulation analysis, we have used 
the 2nd order Gaussian pulse El, 

. . ( 47rf 2 \ /— 27rt 2 \ 

s(t) =A[l- —j exp > ( 2 ) 

The amplitude is adjusted through parameter. A, and pulse 
width is adjusted through parameter, (. 

The received signal is the distorted version of the transmit 
pulse with multipaths. The TOA is defined as the time elapsed 
for the first arrival path to reach the receiver from the trans¬ 
mitter. The received signal can be represented by 

Nf 

Wrx(f) = ^2 d 3 r ^ - i T f ~ C 3 Tc ) + n (*)- (3) 

3 =0 
,- L 

where r(t) = \ jf- X] a l r l(t ~ T l), where Eb, is the captured 
V f i =l 
i=l 

energy and ^ a 2 = 1. The gain and received UWB pulse for 

1=0 

the Z-th tap is given by ai and r;(f). The n{t) is the AWGN 
process with zero mean and double-sided power spectral 
density of Nq/2. With out loss of generality, and for simplicity 


of analysis, we assume Cj = 0 and dj = 1. TOA estimation 
problem is the estimate the first arraival path, ti = T toa , in the 
received signal 0. 

B. System Model 

We define a system model as shown in Fig [Q where 
N receivers each having an energy detector is employed to 
estimate TOA. When the target is moving then the estimated 
TOA is dynamic in nature. In these situations joint fusing and 
tracking using multiple sub-optimal estimates from the sub- 
Nyquist receivers can yield better performance. We employ 
Kalman structure as shown in Fig. |T] to jointly fuse and track 
the TOA. The multiple receivers can be distributed in nature 
and makes independent estimates of the TOA, which are then 
fused using a Kalman filter at the fusion center. All receivers 
should return the TOA for a reference location. Without loss 
of generality, we assume that the all receivers are equidistant 
from the targejj The proposed tracking and fusion strategy 
can also be employed in a single receiver architecture having 
a bank of energy detectors. 

C. UWB Hardware 

There are several low-complexity sub-Nyquist receivers 
based on energy detector available in literature |f25l - l|28l . 
Figure [2] shows a graphical depiction of an in-house developed 
IR-UWB platform for ranging and communication. It uses a 
low cost pulse generator to generate sub-nano second pulses 
using step recovery diode (SRD), as described in (29). The 
characterization and modeling of the UWB platform for a 
distance measurement system can be found in OP) fJD . A 
detailed architectural description and experimental ranging 
results from a prototype of the platform have been published 
in (28) . The power and range of the transceiver can be easily 
traded by controlling the amplitude, duty cycle and number of 
pulses per bit of transmission. Even though proposed methods 
in this paper are not tied to a particular hardware, using 
the above discussed UWB hardware can yield significant 
performance benefits in terms of cost and power since they 
are implemented using low-complexity analog circuits. 

1 If this condition is not met, then each receiver needs to appropriately scale 
the TOA based on the geometry of the reference location and actual location 
of the receiver. 
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(a) Radio module 


(b) Printed circuit board (PCB) 



Fig. 2. Iconic model of the in-house developed impulse radio UWB-platform 
of size 6 cm x 4 cm for ranging and communication working in the 6 GHz 
regime with separate RX and TX antennas from Greenwavescientific. (details 
available at (221) 



T. »(i) 

(a) Matched Filter 



(b) Energy Detector 

Fig. 3. Receiver structure with different detectors. The s(i) and m(n) 
denotes the sampled version of transmit pulse and received signal. 


III. Detectors 

The two commonly employed receiver structures are as 
shown in the Fig. [3] The received signal passes through the 
low noise amplifier (LNA) and band pass filter (BPF) of 
bandwidth B. The output signal of BPF is converted in to 
energy samples, x = (xi, x ?■..., Xn b ), from which TOA is 
estimated. Here, N obs , defines the number of energy samples 
used in the estimation. Based on the energy collection strategy 
different types of detectors exists. In this paper, we consider 
the matched filer which operates at Nyquist rate (1/T S ) and 
energy detectors operating at sub-Nyquist rate ( \/!),). 


A. Matched filter 

In the matched filter, energy is collected by correlating the 
received samples with the transmit pulse shape as shown in 
Fig. [3a] {33). Matched filter can be mathematically expressed 
as 



N p 



= m(n + i)s(i), 

(4) 


i=0 


^-toa 

= argmax(x), 

(5) 

^toa 

= tltoaX's, 

(6) 


where, s(i) = s(iT s ), represent the digitized transmit pulse 
and m(n) = u> a {nT s ), represent the digitized received signal, 
sampled at interval T s . In matched filter, the energy samples, 
x n , are at Nyquist rate, 1/T S . Energy detectors, operating at 
the sub-Nyquist rates are an interesting alternative to digital 
matched filter. 


B. Energy Detector 

The structure for energy detector is as shown in Fig. [3b] 
The structure is amenable for a low-complexity analog imple¬ 
mentation at sub-Nyquist rates m-m- In energy detector, 
the output signal from BPF is converted in to energy samples, 
x = (xi,X 2 ...,XN oh j- TOA is estimated from the energy 
samples using the equation below 


fhm = /(x), 


(7) 


Ttoa — 



( 8 ) 



n 


Fig. 4. Variation of energy samples, x n , verses block index, n. Parameters 
Tf = 200 ns, Tf, = 1 ns, SNR = 0 dB and N 0 = 200 are considered. 


where, /(•), is the estimator function, which estimates the 
block-index/sample-index of the first arriving path. The func¬ 
tion, /(•), is chosen based on the channel model. Equation ©, 
represents f toa as the mid-point of the corresponding estimated 
block, assuming that the true TOA is uniformly distributed 
with in this block. 

Without loss of generality, we consider the sub-Nyquist rate, 
Tb = T c , in detector structure shown in Fig. [3] Therefore, each 
frame consists of, N obs = Tf/T c , blocks, each with energy, 
x n , n £ [1,..., AUs], given by 

N f rU-l)Tf+nT b 

Zn=Y'/ \r(t)\ 2 dt. (9) 

i=0 JU- 1 ) T f + ( n ~ 1 ) T b 

Typical variation of energy samples, x n , verses n is as shown 
in the Fig. Q] In energy detector, the energy samples are at 
sub-Nyquist rate, 1/T b . 


IV. AWGN Channel Analysis 


The best performance in terms of mean-square-error (MSE) 
for an unbiased estimator is given by the CRLB and for a 
TOA estimation problem this is given by EB, G3: 


Gt - 8tt 2 SNR/3 2 ’ 


( 10 ) 


where j3, is the effective signal bandwidth defined by 


njs(/)i 2 4f 


where S(f), is the Fourier transform of the transmit pulse, 

s(t). 
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In AWGN channel model, we consider a single-path model, 
with L = 1, a\ = 1 and n = r t0a in ©. The probability 
density of the x' n , for matched filter and energy detector 
depends on whether the particular energy sample, x n , is 
signal+noise sample or noise-only sample. For matched filter 
it can be shown that the probability density functions under 
hypothesis, Hi = x n is signal+noise sample, and Hq = x n is 
noise-only sample, is given by 


P{x n ) 


J\f(\TEb,a 2 ) under Hi, 
A/”(0, a 2 ) under Hq, 


( 12 ) 


where cr 2 = No/2, is the variance of the noise samples and 
A f denotes the Gaussian distribution. 

To derive the probability distribution function (PDF) of the 
energy samples in energy detectors case, the function /(•) in 
© need to be defined. Under the AWGN assumption, the 
received signal will have a single path, whose delay represents 
the TOA. Thus, optimal TOA estimation strategy here would 
be to pick the energy sample having maximum energy as 
n toa , this selection criteria is called maximal energy selection 
(MES) and is given by 1(211 



Fig. 5. Variation of MSE with SNR for energy detector and matched filter. 
MSE is evaluated by averaging the estimated TOA for 1000 random rt oa 
drawn from U[0,Tf = 200 ns]. Matched filter asymptotically reaches the 
CRLB bound. 


the noise-only energy sample is higher than the signal+noise 
energy sample. That is 


( 21 ) 


fitoa = max(x), (13) 

Aoa “ ^btoa —'j T>, (14) 

The complementary cumulative distribution function (CCDF) 
of energy samples, x n for energy detector, is given by ED, 

ED-Gi 


where x nsn and x nno are energies of signal+noise and noise 
only sample. For AWGN channel, there will be only one 
signal+noise energy sample and all other samples are noise- 
only. The correct selection will happen when x n < x sn , for all 
the energy samples except for the signal+noise energy sample. 
Thus, the probabilty of correct selection, P s , and probalility 
of error selection, P e , is given by f36l 


|Qm '/2 (fs ^r) under Hi Ps 

l ex p(-W?) £ 7 i(tS) underJ? ° Pe 

V x 7 i—0 x 7 


J ( 1 -Q(^-))"' bS \{x n )dx n , 

x n =0 

1 -Ps, (22) 


(15) 

where Qz(cl, b) denotes the Marcum-Q-function with param¬ 
eter Z, P denotes the probability and M’ rj Nf(2BT b + 1), 
denotes the degrees of freedom (DOF). At high SNRs the 
above equation can be approximated to ED 


p(x„) 


M{p Hl ,a 2 Hi ) under Hi, 
N{ph 0 ,Oh 0 ) under H 0 , 


(16) 


The mean and variance of x n under both hypotheses is given 
by 


PH! 

= NfMa 2 + E b , 

(17) 


= 2N{Ma 4 + 4 :<j 2 E 2 , 

08) 

PH 0 

= NfMa 2 , 

09) 

< 

= 2N f Ma 4 , 

(20) 


where M = ‘2!VI), +1, is the degrees of freedom of noise and 
a 2 = Nq/2. 


where p{x n ) is the density function of the energy samples 
for the detectors under hypothesis Hi. Q denotes the CCDF 
function of x n . Here // n() , denotes the mean values of x n 
for matched filter and energy detectors and is given by 0 
and NfMa 2 respectively. a 2 0 , denotes the variance of x n for 
matched filter and energy detectors and is given by cr 2 and 
2NfMa 4 respectively. 

It is not straight forward to arrive at the closed form expres¬ 
sion for energy detector’s P s and P e due to complex distribu¬ 
tion functions (fl5l . There is however a loss of performance 
in energy detectors due to squaring and integration operation 
at sub-Nyquist rates compared to matched filter, operating at 
Nyquist rate. Figure [5j shows the mean square error (MSE), 
performance for matched filter and energy detector. MSE is 
evaluated by averaging the estimated TOA for 1000 random 
TOAs (T toa ) drawn from U[0,Tf = 200 ns . The sampling rate 
(1/T S ) used by matched filter is 1000 GH;^| and sampling rate 
(1/Th) used by energy detector is 1/1000-th that of matched 
filter and is equal to 1 GHz. 


If we choose the MES criteria for TOA estimation for both 
matched filters and energy detectors as given in © - © and 
(fl3l - ( 1 1 41 respectively, then an error can occur, if any one of 


2 Since TOA is continuous in time, very high sampling rate is employed to 
demonstrate that matched filter will indeed reach the CRLB bound without 
much ambiguity due to discretization of time due to sampling. 
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Algorithm 1: Kalman filter for tracking and fusing of TOA 
from multiple sub-Nyquist energy detector estimates 


Input: Prior state, Si n ; t , prior state’s covariance, Mi n i t , and 
covariance of plant noise, Q. 

Output: Tracked and fused TOA estimates, ftoa and variance of 
estimate, <f t o a 

«•[!]+- 

M[l] -r- ivj.j n i t 


M m 


> Initial state 
Initial states Covariance 

Initial TOA estimate 
Initial MSE estimate 

> Pridicted state 

> Predicted MSE 


Ttoa[l] +" s init[l] > 

<Ttoa[l] M; n i t [l](l, 1) > 

for i +- 1 to n do 

Sprd As[i] 

M prd +- AM[t]A' + s 

> Compute Kalman gain 
K = (M prd H T )(C + HMprfH 1 )- 1 

> Update state from observations {y[n]) 
s[i + 1] = s prd + K(y[i] - Hs prd ) 

Ttoa[* + 1] = s[i + 1] (1) > Save fused TOA (ftoa) 

M[i + 1] = (I 2 - KH)M prd > Update MSE 


^toa \} + 1] — M[i + 1] (1, 1) 


> Save MSE 


13 return f toa ,S>1 


At high SNRs, to achieve the performance of the matched 
filter using multiple energy detectors, we need the covariance 
of the signal+noise energy sample (under hypothesis Hi) to 
be same. We can approximate the probability distribution 
of energy samples for energy detectors in to a Gaussian 
distribution at high SNR as shown in (fl6t . If there are K 
receivers having energy detectors to independently estimate 
TOA, then signal+noise block in each detector will have a 
probability distribution given by 


G ~ A/ - , cr^) ,i £ [!,■■• , K]. (23) 


Since e,.s are independent and identically distributed, then 
the best estimate for the signal+noise block using K energy 
detectors is given by ED , 


e 


! 

Ned 


J2 ek - 


(24) 


The probability distribution function of e is given bj0 


' N 


Mjtd 


U H i 

K 


(25) 


Let er^p = a 2 and a^ D = crjj denote the variance of the x n 
under Hi for matched filter and energy detector. From ( 1 1 21 
and (l25l > to achieve the performance of the matched filter using 
multiple energy detectors, we need 


J MF 


A+d 


'EDi 


(26) 


where, Ned, is the number of energy detectors needed to 
have the same performance as that of matched filter. For 
normalized energy per bit (Ef, = 1), and using the variance 
for signal+noise sample from (fl2l >. ( 1 1 Gb and ( | I S[ i we get 


Aed 

= lint 

(27) 


SNR—>-oo ^mf 

Aed 

(2 Mcr 4 + 4cr 2 ) 

= lim 

(28) 


a 2 —X) (J 2 



= 4. 

(29) 


Thus, from ( l29] i. asymptotically, with the increase of SNR 
(Eh/No), the number of energy detectors needed to achieve 
the same performance as digital matched filter is equal to 4. 
We will show in simulation that this phenomenon is indeed 
true in the later section. 


V. Multipath Channel Analysis 

Many UWB ranging applications have channel response 
with several multipath components, i.e. the received pulse in 
© has (or, ct 2 ■ ■ ■ocl',t~i,T 2 ■ ■ ■ ,tl), where L is the number 
of multipaths. The TOA estimation problem is to identify 
the leading edge (first arriving path, n). In multipath UWB 
channels, the matched filtering performance is not optimal 
since the shape of the pulse is lost in the channel due to 
the frequency selective fading. Also, the magnitude of the 

3 See (38) to obtain the probability distribution function of a function of 
Gaussian random variable. 



Fig. 6. Fusion of multiple Detectors using Kalman Formulation. 
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Fig. 7. Variation of MSE with SNR for energy detector (with selection criteria 
MES and WMESS) and matched filter. MSE is evaluated by averaging the 
estimated TOA for 1000 random rtoa drawn from IA[ 0,100 ns]. 


Fig. 8. Averaged Channel energy profile for 802.15.4a CM1 channel model. 
100 different channel realizations are averaged. 


energy sample containing first arriving path may be smaller 
than the peak energy sample, therefore, using MES criteria 
for TOA estimation, as represented in ( IT3l ) and (J5) does not 
always yeild true TOA (39) . Figure [7J shows the performance 
in terms of MSE for matched filter and energy detector for 
802.15.4a, residential LOS channel (CM1 model) (39). MSE 
is evaluated by averaging the estimated TOA for 1000 random 
T t0 a drawn from U[ 0,100 ns]. The sampling rate (1 /T s ) used 
by matched filter is 8 GHz and sampling rate (1/T&) used by 
energy detector is 1 /8-th that of matched filter and is equal to 
1 GHz. 

The performance of the energy detectors can be improved 
by considering the a-priori information such as power delay 
profile (PDP). To accomplish this, we can repose the TOA 
estimation problem as a multiple hypothesis testing problem. 
The energy samples vector, x, is a 1 x N 0bs vector, out 
of which, N e blocks are having multipath signals, N e T b , 
represents the excess delay of the channel. If we define Hk, 
as the hypothesis that the k-th sample denotes the first path 
arrival, then Hk = x, with. 


x(n) = < 


rnT b 

/ n 2 (t).dt , 

J(n-l)T b 

for n = 1,..., k — 1 

r n T h 

/ |wrx(f) + n(t)\ 2 .dt, 

J(n-l)T b 

for n = k,..., k + N e — 1 

r n T b 

/ n 2 (t).dt , 

J(n-l)T b 

for n = k + N e ,..., N ohs . 


(30) 


Here, x(n), denotes the n-th element of x and H nt , is the 
true hypothesis ED. If the PDP for the channel is available, 
then we can modify the hypothesis test from MES to weighted 
maximum energy sum selection (WMESS), and is given by 


fitoa = arg max \x(k : k + N e )£} , (31) 

fce[i,...,JV obs ] 

where, x(k : k+N e ) is a 1 x N e vector having the N e elements 


from x starting from k and £ is the N e x 1 vector denot¬ 
ing the a-priori channel energy information which is similar 
to correlating received energy samples with PDP. Averaged 
channel energy profile, £, for IEEE 802.15.4a residential LOS 
channel model is as shown in Fig. [8] The performance of 
energy detector with WMESS algorithm is as shown in Fig. [7] 
Even though, the WMESS estimation algorithm require more 
computation, its performance is better than matched filter at 
low SNRs for IEEE 802.15.4a CM1 multipath channel model. 
However, at high SNR there is a loss of performance compared 
to matched filter. 

The performance degradation at high SNR with WMESS 
criteria for energy detector can become a problem when the 
TOA of the estimated target is dynamic in nature. In the next 
section, we will discuss a mechanism to fuse the estimates 
from multiple receivers to achieve improved performance for 
a dynamic TOA system. 


VI. Multi detector Fusion Using Kalman Filter 


When the estimated TOA of the target is dynamic in nature, 
then joint fusing and tracking using multiple independent 
TOA estimations from the energy detectors can yield better 
performance. We employ Kalman structure as shown in Fig. |T] 
to accomplish this. 

The state of the system and its dynamic nature is represented 
by the state equation in the Kalman filter formulation, and is 
given by 

s[n] = As[n — 1] + u[n], (32) 


where s[n]. A, and u[n] is given by 


s n = 


AW = 


T[n\ 

Vi oa [n] 

1 1 ' 
0 v 


(33) 

(34) 


u[n] 



(35) 


where u v [n\ is the noise with variance a 2 and v < 1 is a 
constant. The r[n] is the TOA at n-th time interval and u toa is 
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an AR(1) processes with constant v, representing the dynamic 
nature of TOA. The plant noise u[n] ~ M (0, Q), where Q is 
given by 


Q = 


The measurement matrix, y = 


u p\ 



[^toa ? 

T 2 

1 toa> 


'N' 


(36) 


, is the in¬ 


dependent TOA estimates from N receivers. The measurement 
equation is given by 


where 


y[n] = H[n]s[n] + w[n], 


H[n] 


1 0 
1 0 


1 0 


(37) 


(38) 



Fig. 10. Steady state MSE verses time, for multiple energy detectors and 
digital matched filter. 


w[n] ~ J\f( 0, diag (ct 2 ,... &%)). N denotes the number of de¬ 
tectors and diag(-) indicates the diagonal matrix with diagonal 
elements represented inside the brackets. 

The Kalman filter based fusion of dynamic TOA from 
multiple energy detector estimates, y[n], is as shown in 
Algorithm!]] At each iteration, the Kalman filter yields fused 
TOA estimate, r toa [n], and its variance, a t 2 a [n], from multiple 
energy detector estimates, y[n[. In our analysis, we have used 
same kind of detectors and also the reliability or importance 
of each of the detectors are assumed same. By appropriately 
selecting the measurement matrix H and covariance of the 
measurements, we can extend the Kalman filter formulation to 
accommodate different types of detectors with varying char¬ 
acteristics. In the next section, we will assess the performance 
of the proposed fusion and tracking method. 



Fig. 11. Number of sub-Nyquist energy detectors need to match the 
performance of matched filter operating at Nyquist rate. 


VII. Simulation Study 
A. AWGN Channel 

AWGN channel is simulated with F s = 8 GHz with 
Tf = 200 ns, N s = 1 and B = 4 GHz. We set X& = 1 ns, 
thus resulting sampling rate of the energy detector is 1 /8-th 
the Nyquist rate. Consider a linear variation of TOA (r toa ) 
observed by a set of four independent energy detectors (refer 
Fig-03- The observed TOA of the four detectors are as shown 
in Fig. [6a] at SNR = 12 dB. The tracked TOA using the 
Kalman formulation is shown in the Fig. [6b] We consider prior 
state, covariance of prior state and plant noise variance equal 
to [20,1], O.OII 2 and cr 2 = 0.0001 respectively. 

The mean square error (MSE) variation with time for the 
detection schemes with multiple energy detectors at various 
SNRs ( Eb/No ) are illustrated in the Fig. [9] Notice that the 
steady state variance decreases as the number of energy 
detector estimates increases. Also, from Fig. [9] at high SNRs, 
only few energy detectors (operating at 1 /8-th rate of matched 
filter) are sufficient to achieve same performance as digital 
matched filter. 

Figure [TOj shows the variation of steady state MSE with 
number of energy detectors (operating at F s / 8). Notice that 
steady state MSE of fusing 4 independent energy detector 
estimates reaches that of digital matched filter estimate. 


Figure [TT] shows the number of energy detectors operating 
in sub-Nyquist (F s / 8) rate needed to match the performance 
of matched filter operating at Nyquist rate ( F s ). At lower 
SNRs, more number of energy detectors, (IV ED ), are required 
to meet the matched filter performance. The number of en¬ 
ergy detectors, Ned, reduces with the increase of SNR and 
asymptotically approaches to 4, confirming with the analytical 
derivation of the previous section. 

B. Multipath Channel 

IEEE 801.15.4a CM1 channel model is simulated with 
F s = 8 GHz with T f = 200 ns, N s = 1 and B = 4 GHz. 
We set T;, = 1 ns, thus resulting sampling rate of the 
energy detector is 1 /8-th the Nyquist rate. We use WMESS 
algorithm discussed in the previous section to arrive at TOA 
estimates from the energy samples. We employ similar fusion 
and tracking strategy using the Kalman filter described in the 
previous section for the dynamic TOA model. The prior state 
information, covariance of prior state and other initialization 
parameters of the Kalman filter are kept same as in the 
previous section. Figure [T2] shows the variation of MSE 
with time at SNR of 20 dB for different number of energy 
detectors. As expected the steady state variance decrease with 
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- Ned = 1 

- Ned ~ 4 

- N ed = 7 

— Ned = 10 

Mached Filter at Nyquist Rate 




(c) MSE verses time(SNR = E^/Nq = 10 dB) (d) MSE verses time (SNR = E^/Nq = 20 dB) 


Fig. 9. Variation of MSE with time for digital matched filter and fusion of multiple sub-Nyquist energy detectors using Kalman filter at different SNRs 


the increase in the number of detectors. Figure [13] shows the 
variation of steady state MSE with number of sub-Nyquist 
energy detectors, as discussed earlier a single energy detector 
with WMESS selection criteria outperforms the matched filter 
at low SNRs, however, at high SNRs more energy detectors 
are need to match the matched filter performance. Figure 
\n\ also illustrates that fusing two energy detector estimates 
(Ned = 2) can match the matched filer performance across 
the SNR region. 

VIII. Conclusion and discussion 

Estimating TOA with good accuracy is very important for 
localization and several other applications. UWB with its 
large bandwidth can yield high precision ranging as evident 
from CRLB Equation For AWGN channel, matched 
filter based estimation yields optimum performance at high 
SNRs. However, matched filter requires Nyquist-rate sampling 
and is inefficient in terms of cost and power consumption. 
Energy detectors operating at sub-Nyquist rate is an interesting 
alternative as it can be designed using cost effective analog 
circuits and is power efficient, however, it lacks the precision in 
range measurements. In this paper, we showed that for AWGN 



t [s] 


Fig. 12. Variation of MSE for IEEE 802.14.4a CM1 channel model, with time 
for digital matched filter and fusion of multiple energy detectors (WMESS 
criteria is employed). 


channel, we can achieve the performance of a matched filter, 
by using a multiple sub-Nyquist energy detector estimates. 
We showed that number of energy detectors needed to meet 
the matched filters performance is high at low SNRs and 































































9 



Fig. 13. Steady state MSE verses time for IEEE 802.14.4a CM1 channel 
model, for multiple energy detectors (WMESS criteria is employed) and 
digital matched filter. 

reduces as SNR increases, and finally converges to 4 as SNR 
increases asymptotically. This is analytically derived in < [29t . 
and confirmed through simulations (refer Fig. [TO] and Fig.ITTli. 
A Kalman filter with suitable choice of state-equation and 
measurement equations is designed to perform the dual task 
of tracking the TOA as well as fusing the multiple energy 
detector outputs. Filter equations are shown in Algorithm Q] 
and performance in terms of MSE are demonstrated in Fig. [9] 
and Fig. [TO] Result indicate that for AWGN channel, the steady 
state variance drops with the increase of number of detectors 
and requires 4 energy detectors to have the same performance 
as digital matched filter. 

Proposed tracking and fusion strategy of energy detectors 
is analyzed in simulation for IEEE 802.15.4a CM1 multipath 
channel (residential LOS model) and results are demonstrated 
in Fig. [7] Fig. [12] and Fig. [13] Results indicate that the two 
sub-Nyquist sampled energy detector estimates with WMESS 
criteria can outperform digital matched filter at all SNR 
regions. 

We plan to further this work, by assessing the performance 
of the proposed fusion and tracking using multiple energy 
detectors using our in-house UWB radio hardware discussed 
in Section HU 
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